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ABSTRACT 

This  paper  examines  the  standard  methods  for  analog  differentiation. 
It  concludes  that  the  compromised  derivative  circuit  with  paralled  cap- 
acitance and  resistance  in  the  feedback  loop  and  added  resistance  in  the 
input  leg  is  the  most  satisfactory  normal  circuit  provided  that  the  poles 
generated  by  the  components  are  suitably  located.   It  shows  furthermore 
that  multipole  derivative  circuits  can  be  devised  from  mathematical 
approximations  of  the  type 

.m  (m) 

Such  circuits  have  greater  descrimination  against  frequencies  above  a 
specific  level  and  give  more  accurate  derivate  for  frequencies  which 
are  low  relative  to  this  cut-off  point. 
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I .   INTRODUCTION 

The  importance  of  quantitative  description  has  been  well  stated  by 
Lord  Kelvin. 

"I  often  say  that  when  you  can  measure  what  your  are 
speaking  about ^  and  express  it  in  numbers,  you  know 
something  about  It;  but  when  you  cannot  measure  it, 
when  you  cannot  express  it  in  numbers,  your  knowledge 
is  of  a  meagre  and  unsatisfactory  kind;  it  may  be  the 
beginning. of  knowledge,  but  you  have  scarcely,  in 
your  thoughts,  advanced  to  the  stage  of  science, 
whatever  the  matter  may  be . " 

Today  in  experimental  engineering,  in  system  and  process  control, 
we  are  concerned  with  the  measurement  and  use  of  electrical  signals 
derived  from  a  wide  range  of  varying  physical  states.   Transducers  exist 
which  can  convert  the  variation  of  measurable  quantitites  into  such 
signals.   However,  at  the  present  time,  with  the  exception  of  velocity 
type  pick  ups,  there  are  no  transducers  which  give  outputs  proportional 
to  the  time  rate  of  change  of  the  transduced  property.   Likewise  in 
many  systems  used  for  process  control  a  derivative  feed  back  loop  is 
preferable  to  an  integral  feed  forward.   The  factors  of  primary  impor- 
tance in  connection  with  rate  mode  are  that  by  opposing  all  change,  rate 
mode  has  a  great  stabilizing  effect  on  control. 

The  need  for  a  clear  understanding  of  electrical  analog  derivative 
circuits  is  apparent.   On  the  assumption  that  the  signal  input  is  clean, 
standard  text  books  on  analog  methodology  give  circuits  which  will  pro- 
duce derivative  outputs.   However,  the  world  of  reality,  as  we  know  it, 
is  always  noisy,  and  consequently  this  theoretical  ideal  input  is  never 
achieved.   Thus,  the  derivative  circuit  responding  to  the  noisy  input 
differentiates  both  the  noise  and  the  pertinent  signal.   Since  in  many 
natural  phenomena  and  control  processes  the  variable  changes  are 
restricted  to  relatively  low  frequency  and  noise  is  primarily  generated 
or  picked  up  at  high  frequency  the  signal  to  noise  ratio  is  worse  after 
differentiation.   This  is  due  to  the  fact  that  a  derivative  device  is  by 


Superscript  numerals  refer  to  the  Bibliography  at  the  end  of  this 
report . 
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definition  a  rate  of  change  responding  device.   Standard  authorities 
make  it  clear  that  the  very  nature  of  the  response  and  of  the  electronic 
circuitry  makes  this  unsatisfactory  situation  one  which  the  process 
control  and  instrumentation  engineers  must  learn  to  live  with,  to  avoid 
or  to  compromise.   Engineering  is  of  course  the  art  of  compromise  —  but 
for  such  to  be  effective  its  implication  must  be  understood.   The  purpose 
of  this  paper  is  to  attempt  to  demonstrate  how  analog  derivative  circuits 
can  be  adapted  to  specific  performance  requirements. 


II.   STAMDARD  DERIVATIVE  CIRCUITS 


The  circuits  given  in  Figs  la-c  are  those  generally  quoted  "by 
authoritative  texts  as  the  most  suitable  for  analog  differentiation.  '    ' 
The  circuit  shown  in  Fig.  la  is  a  true  differentiator.   This  is  readily 
see  from  the  system  analysis  given  below. 

e^{s)   =   -SRCe.^(s)  (l) 

which  may  be  written  in  the  time  domain 

e„(t)  =  -RC  •  |-  e.   (t)  (2) 

0  dt   m 

This  latter  equation  is  of  course  a  special  case  of  the  more  general 
Laplace  transform  theorem. 


^(g) .  .>.,  - 


vdt"/  ,^,   dt^"^ 

k=l 

The  validity  of  Eq.  (2)  is  not  dependent  upon  the  form  of  the  input  signal, 
it  merely  requires  that  n  =  1.   Thus,  the  circuit  will  not  discriminate 
against  any  frequency,  but  will  give  the  derivative  of  the  input  which  it 
receives.   Let  us  consider  then  its  behavior  when  the  input  is  defined  by 

e.   =  A  sin  wt  (k) 

m 

e    =  -coRCA  cos  wt  (5) 

out 

¥e  see  straightaway  that  the  sinusoidal  signal  is  phase  shifted   7r/2 
and  increased  in  amplitude  by  the  factor  wRC.   Consider 

e.   =  f(t)  +  A  sin  wt  (6) 


and 


e    =  f '(t)  +  ojA  cos  wt  (7) 


It  is  clear  that  unless  f(t)   is  some  simple  function  of  wt  the 
relative  amplitudes  of  the  two  components  are  not  the  same  in  the  initial 
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FIG.  I.      DERIVATIVE  CIRCUITS 
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function  and  the  derivative.   However^  in  the  particular  case  for  which 
f(t)  =  e     we  meet  the  special  condition  that 

-e  ^  =  (jof(t)  +  wA  cos  cot  (8) 

out 

In  analog  control  loops  the  pertinent  signals  are  generally  restricted 
to  the  low-frequency  range  but  this  essential  data  is  almost  always  mixed 
with  "noise".   Noise  is  essential  high  frequency.   Thus  the  output  of 
this  true  differentiation  circuit  will  be  the  desired  derivative  with  the 
noise  derivative  superposed  and  by  virtue  of  the  fact  that  the  special 
case  of  Eq.  (8)  is  not  usually  met;  the  signal  to  noise  ratio  will  be 
seriously  reduced.   Indeed,  in  many  applications  a  true  differentiator 
becomes  a  noise  amplifier,  e.g.,  when  f(t)   is  zero. 

The  circuit  of  Fig.  lb  is  designed  to  solve  the  implicit  differ- 
ential equation 

(l-a)i..=i  (9) 

and  so  to  obtain  an  approximation  to  the  derivative.   For  the  special 
case  when  (a)  is  unity  the  actual  derivative  results.   Clearly  in  this 
instance  there  can  be  no  advantage  in  this  circuit  since  it  is  more 
complex  than  that  of  Fig.  la  and  gives  an  identical  output.   It  is  the 
case  a  ^  0,  which  is  of  interest.   The  above  equation  is  implemented  on 
a  computer  by  first  integrating  both  sides  of  Eq.  (9)^  resulting  in 

t 
(1  -  a)z  +  /   zdt  =  X  (10) 

0 

or 

t 

z  =  az  -  /   zdt  +  X  (ll) 

0 

The  transfer  function  of  this  circuit  is 

z         S 


X   "  (1  -  a)S  +  1 
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(12) 


or  -with   T  =  (l  -  a)  we  may  write 

f  =  -  (St  +  l)  ^^^^ 

It  is  seen  from  Eq.  ( 12 )  or  (l3)  that  when  a  =  1,  the  true  deriva- 
tive results,  as  we  previously  noted.   Figure  2.,    shows  the  frequency 
response  In  this  case,  (llnel).   For  the  general  case  a  7^  1.   The  fre- 
quency response  is  limited  for  circular  frequencies  greater  than  w 
defined  by 

1   T   (1  -  a)  ^      ' 

This  is  illustrated  in  Fig.  2  (line  2).   Thus,  we  readily  appreciate  that 
the  introduction  of  the  single  pole  causes  input  signals  of  frequency 
higher  than  co,   to  pass  with  a  constant  amplification.   Hence,  we  may 
achieve  a  considerable  improvement  in  signal  noise  ratio  in  a  particular 
circumstance  by  the  correct  positioning  of  the  pole,  i.eo,  by  an  appropriate 
choice  of  the  parameter  (a). 

The  transfer  function  for  the  circuit  given  in  Fig.  Ic  is 

6^(5)  R^C^S 

e.^(s)  =   ■  (R^C^S  +  IJlR^CgS  +  1)  ^ ^^ ^ 

For  the  special  case  in  which  R   and  Q,        are  zero  and  RpC-,  ^  0  this 
equation  is  identical  to  Eq.  (l)  and  the  discussion  of  circuit  l(a)  is 
relevant.   In  the  special  case  when  either  R,   or  C^   is  zero,  the 
transfer  function  is  similar  to  that  for  l(b)  and  the  behavior  then 
follows  the  preceeding  discussion.   When  these  particular  conditions  are 
not  met  the  transfer  function  contains  the  zero  necessary  to  generate  the 
derivative  together  with  two  arbitrarily  located  poles.   For  convenience, 
at  this  stage  of  our  discussion,  we  shall  put  R-,C-,  =  Rp'^o  -   '^ '      I^  this 
case  Eq.  (l5)  niay  be  re-written  as 


e  (S)      RC  S 
°  2  ^  (16) 


e.  (S)     i^  ^2 

irt  '  (St3_  +  1) 
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The  high  and  low  frequency  asymptotes  of  this  transfer  function  are  given 

in  Fig.  2  (line  3).   The  circuit  has  the  advantage  of  passing,  with  no 

2 
amplification,  frequencies  of  go  ,  and  attenuating  all  higher  frequencies. 

If,  for  design  purposes,  it  is  desired  to  operate  upon  signal  frequencies 

up  to  w   with  a  maximiom  error  of   1  per  cent,  the  time  constant  values 

should  be  chosen  such  that 

Since  the  error  decreases  approximately  as  the  square  of  the  frequency 
difference  from  co  ,  at  a  frequency  of   n/ 10  oo    the  error  will  be 
approximately  0.1  per  cent. 

The  frequency  response  diagram  given  in  Fig.  2  shows  clearly  that 
it  is  possible  by  the  use  of  multiple  poles  to  cause  a  derivative  type 
circuit  to  discriminate  against  certain  frequencies  and  yet  maintain 
reasonable  accuracy  in  a  restricted  bandwidth.   Thus,  the  multiple  pole 
system  has  distinct  advantages  in  control  applications.   Generally 
speaking,  a  control  signal  consists  of  either  a  dc  or  low  frequency 
signal  and  a  superposed  higher  frequency  noise.   In  most  laboratories, 
the  critical  noise  is  60cps.   Thus,  in  a  normal  laboratory  control 
application  the  circuit  l(c)  finds  the  most  common  acceptance.   When  it 
is  used  it  is  customary  to  require  that  the  60  cycle  noise  passes  without 
amplification.   The  location  of  the  poles  is  then  determined  in  accor- 
dance with  the  following  proceedure. 

At  00  =  377  rad/sec,  the  gain  is  zero  db.   Hence 

=  0  (18) 


(ST  +  if 


or 


db 


(jO 

=  1  (19) 


y  [1  -  ((.t)2]2  +  (2c.t)2 


therefore 


[1  -  (a)T)2]2  +  (2cot)2  =  ^2l  (2o) 

1  -  2(aJT)2  +  (oot)^  +  4(u)t)2  =  u)2  (2l) 

[1  +  (a)T)2f  =  0)2  (22) 


i.e 


1  +  (a)T)2  =  to  (23) 

Hence 

T  =  —  =  .0515  sec  (24) 


Thus 

CO  =  —  =  19.4  rad/sec  (25) 

It  can  "be  shown  that  if  the  corner  frequency  is  fixed  at   19-4  rad/sec 
a  1  per  cent  error  exists  at  1.9^  rad/sec  or  0.31  cps.   Thus  the  circuit 
discriminates  against  the  higher  frequencies  at  the  expense  of  accuracy 
in  even  the  very  low  frequency  spectrum.   From  Fig.  2  it  is  clear  that 
with  a  two  pole  approximate  derivative  circuit  it  is  preferable  to  make 
the  poles  coincident. 
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III.   CIRCUITS  BASED  ON  MATHEMATICAL  APPROXIMATION  TO  THE  DERIVATIVE 


3.1.   Feasabllity  of  Approach 

Because  of  the  practical  problems  experienced  in  derivative  circuits 
it  was  decided  that  it  might  be  profitable  to  investigate  the  feasibility 
of  employing  mathematical  approximations  in  the  construction  of  such 
circuits. 

The  derivative  of  f(t)  with  respect  to  t   is  defined  as 

f'(t)  =  ^^"^  ^^^      ^  (26) 

At^O      ^^ 

where  f   is  the  value  of  f(t)   at  time   (t)   and  f   is  value  at 
time  (t  -  At).   The  ratio   (f   -  f  )/At  approximates  to  the  derivative 
when  At   is  small,  the  smaller  At  the  better  the  approximation. 
Physically,  this  fundamental  definition  can  be  interpreted  in  the 
following  way: 

The  time  derivative  of  a  varying  voltage  is  equal  to  a  constant 
multiplied  by  the  change  in  signal  amplitude  which  occurs  in  a  small 
but  fixed  interval  of  time. 

We  examine  then  the  nature  of  this  approximation  to  a  derivative 
when  the  system  is  applied  to  a  signal  which  consists  of  two  components  — 
a  smooth  signal  and  a  superposed  AC  signal.   If  the  time  interval  between 
the  sample  f   and  the  sample  f_   is  correctly  chosen,  the  AC  components 
will  be  in  phase,  equal  in  magnitude  and  of  the  same  polarity.   Thus, 
when  the  signal  samples  are  differenced  the  AC  components  will  cancel 
each  other  out  and   (f   ~  ^q)  will  be  approximately  proportional  to  the 
derivative  of  the  smooth  component.   We  have  implied  in  this  approach 
that  the  sampling  period  must  be  a  function  of  the  frequency  of  the 
AC  signal.   The  reader  will  notice  that  there  is  a  direct  comparison 
between  this  approach  and  Lanczo's    differentiation  technique  for  a 
Fourier  Series.   It  is  well  recognized  that  if  we  differentiate  a 
Fourier  Series  we  usually  get  what  is  not  wanted  in  the  physical  problem, 
namely,  a  large  high  frequency  oscillation.   This,  of  course,  is  precisely 
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what  happens  when  we  use  a  classic  approach  to  obtain  the  derivative  of 
a  noisy  signal.   Lanczo  introduced  a  practical  approach  to  such  a 
problem  when  the  suggested 

—  (27) 

27r  ^  '^ 

as  an  estimate  of  the  derivative  of  f(t)   at  the  point  t   instead  of 
the  limit  process  which  has  dubious  physical  significance  when  we 
differentiate  a  Fourier  Series  approximation. 

The  question  we  must  now  answer  is  how  can  this  simple  approach  be 
embodied  in  a  computing  circuit.   There  appears  to  be  two  distinct  possi- 
bilities.  Firsts  we  can  sample  the  signal  at  discrete  intervals  of  time 
and  difference  the  samples.   Alternatively,  we  can  bifurcate  the  signal, 
delay  one  component  relative  to  the  other  and  then  recombine  the  two 
parts  through  a  summing  network.   The  output  of  this  network  will  be 
proportional  to  the  derivative. 

It  seems  apparent  that  the  first  approach  is  more  suited  to  a  digital 
system  and  the  second  to  an  analog  method.   It  is,  therefore,  to  the 
latter  that  we  shall  direct  our  attention.   We  begin  by  considering  how 
this  approximation  fits  into  an  electronic  network  and  the  first  point 
which  we  clearly  observe  is  that   l/At   is  the  overall  circuit  gain. 
Thus,  the  system  does  not  allow  complete  freedom  in  determining  At  but 
rather  At  must  be  chosen  to  be  compatible  with  the  characteristics  of 
the  amplifiers  which  are  used.   We  must,  of  course,  also  consider  the 
response  characteristics  we  desire  when  assessing  the  circuit  values. 
A  cursory  examination  would  indicate  that  only  one  time  delay  network 
is  required  to  generate  the  potential  difference  which  is  proportional 
to  the  derivative.   However,  if  the  best  discrimination  against  noise  is 
required,  this  is  not  the  case. 

3.2.   Two  Pole,  Two  Leg  Derivative  Circuit 

The  analog  circuit  used  is  shown  in  Fig.  3.   The  transfer  function 
of  the  circuit  is,  therefore 
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FIG.  3.     TWO  POLE,  TWO  LEG  DERIVATIVE  CIRCUIT 
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v^ 


'^1   r^  +  R 
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r^  +  R^ 


m)  ■  I 


R, 


a,  = 


1   r^  .  R^ 


R. 


f2 


R, 


•N 


«a 


r  +  R 
2    2 


J 


(28) 


°2"g"2       „    ffg    "f 
'2  -  r^  +  Rg   '    2  -  R  •    r^  +  Rg 


The  above  equation  can  be  written  as 


e.  (s)  "  ~^\^^Z   ~   ^1^    (St,  +  1)(St^  +  1) 


(29) 


if  the  parameters  are  chosen  such  that 


and 


^1  =  ^=^2 


T^^  T. 


(30) 
(31) 


With  the  exception  of  the  gain  factor  OiAi     -   f  ),  which  must  be 
compensated  for  in  the  summing  amplifier^  the  analysis  of  this  circuit 
is  similar  to  that  for  two  pole  approximation  discussed  previously.   The 
difference  is  due  to  the  fact  that  the  two  poles  must  be  separated  and 
hence  the  desired  values  of  the  corner  frequencies  can  best  be  obtained 
by  graphical  methods  using  the  asympotes  of  the  frequency  response  gain 
curve . 


Accuracy  dictates  that  the  poles  should  be  closely  spaced,  but,  the 
■  ''•-I  )  gain  factor,  dictates  their  values.   This 
that  mathematical  approximations  may  have  advantages 


(x   -  T  )  gain  factor,  dictates  their  values.   This  approach  demonstrates 


-  13 


3.3.   The  General  Case 

In  a  paper  published  in  the  Math  Gazette  ,      W.  G.  Bickley  has  given 
a  set  of  tables  for  use  in  numerical  computation  of  derivatives.   In 
these  tables  coefficients  are  given  for  insertion  into  the  general 
formula 


i-  =  fT  (Vo*  Vi"  ••■  Vp'  *^  f'^^' 


m:    f;  ^  u"u   1"!      p-'p 

In  these  formulae,  m  =  1,  2.,    3,  etc.,  and  reference  to  the  order  of 
the  derivative  of  y  with  respect  to  x  at  the  point  x  =  ra  where 
r  =  0,  1,  2,  ....  p  . 

It  is  significant  to  note  that  in  these  general  formulae  the  con- 
dition is  always  met. 


n=p 


A  =  0  (33) 

n 


n=0 


This  general  formula  indicates  that  if  we  sample  a  signal  at  discrete 
times  At,  2At;  etc.,  amplify  these  samples  by  an  amount  proportional  to 
Bickley 's  coefficients  and  algebraically  sum,  the  sum  so  derived  is  pro- 
portional to  the  n  "  derivative.   The  number  of  time  dealys  (At's)  used 
in  the  process  determines  the  resulting  number  of  poles  generated  in  the 
overall  transfer  function  of  the  circuit. 

3.4.   The  Three  Pole,  Three  Leg  Derivative  Circuit. 

From  Eq.  (32),  the  formula  for  the  three  point  approximation  to  a 
derivative  is 

Figure  \   gives  the  analog  implementation  of  the  above  equation. 
Setting 
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FIG.  4.      THREE  POLE  ,  THREE  LEG  DERIVATIVE  CIRCUIT 
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°l'-l«l 
^l'\ 


R. 


a    = 


1   r^.R^ 


T^  = 


"3  "  r,  +  R. 


the  transfer  function  of  the  circuit  becomes 


e,(s) 


^A 


^2^2 


^3^3 


St  +1   St  +  1   St^  +  1 


(35) 


Based  on  the  selection  of  the  time  constant  such  that   t  =  nT  as 

n 

suggested  by  Bickley's  constant  interval,  the  transfer  function  can  be 
written  as 


e„(s) 

7 — r  =  2a-a^T^ 

e  .  (  s  j     111 
m 


(St^  +  IJiST^  +  1) 


(36) 


provided  the  circuit  parameter  are  chosen  such  that  the  following  con- 
ditions are  satisfied: 


^1^1  -  ^2^2  ^  ^3^3  =  ° 


6a  a  -  3a  a  +  2a  a  =  0 


5a^a^  -  4a20;2  +  Sa^a^  ^  0 


(37) 
(38) 

(39) 


The  ratio  that  satisfies  the  above  three  equations  was  found  to  be 
1,  -k,    3.   This  ratio  is  identical  to  the  coefficients  listed  by  Bickley 
for  his  three  point  differentiation  formula. 

Thus,  a  three-pole  approximation  has  been  generated  using  the 

mathematical  finite  difference  method.   The  value  of  three  poles  over 

two  is  obvious,  they  permit  greater  accuracy  in  the  lower  frequency 

range.   The  net  effect  of  three  poles  can  be  a   -12  db/octave  slope 
through  the  higher  frequency  chosen  to  have  0  db  gain. 
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Figure  5  shows  the  theoretical  frequency  response  of  the  three-pole 
approximation  derivative  circuit  given  in  Fig.  h.      The  time  constants  used 
were 

T   =  .011,   T   =  .022,   and   t   -  .033   . 

The  estimated  response  of  the  double-pole  circuit  with  equal  time  con- 
stants is  included  for  comparison.   The  curves  illustrate  clearly  that 
the  three-pole  circuit  has  considerably  more  accuracy  in  the  lower  fre- 
quency region,  although  both  circuits  pass  60  cps  without  gain. 

3.5-   The  Four-Pole,  Four-Leg  Derivative  Circuit 

The  improvement  of  the  low-frequency  response,  resulting  from  the 
inclusion  of  an  additional  pole  leads  one  to  consider  the  practicality 
of  a  four-pole  derivative  network.   Following  the  procedure  used  in  the 
previous  sections,  the  desired  derivative  is  expressed  in  the  form 

y'  =  K(A^y^  +  A^y^  +   A^y^  +  A^y^)  (ko) 

Figure  6  shows  the  analog  implementation  of  this  equation.   Setting 
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the  transfer  function  of  the  above  circuit  becomes 
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Again  chosing  the  time  constants  such  that   t  =  nT(n  =  1,2, 3,^),  the 
transfer  function  becomes 


m 


(St^  +  DlST^  +  DCSt^  +  1)(St^  +  1) 


(^2) 


if  the  circuit  parameters  are  chosen  such  that 


12a^a^  -  ^a^OL^  +   ^a^a^  -  3a^a^  =  0 


26a  a  -  19a  a  +  l4a  a  -  11a,  a  =  0 


^A  "  ^2^2  ^  ^3^3  "  %\   "  ° 


9a-^Q!^  -  Sa^a^  +  Ta^a^  -  ^W  f   0 


(^3) 
(^6) 


The  solution  of  the  above  equations  yields  the  following  relation- 
ship between  the  parameters 

{hi) 


^2^2  =  12a^a^ 


^3^3  =  2Tai«i 


%\  =   i6  a^a^ 


(^8) 
(^9) 


At  this  point  it  must  be  noted  that  the  ratio  among  the  parameters  as 
given  by  Eqs .  (^7)-(^9)^  is  not  the  same  as  the  coefficients  listed  by 
Bickley  for  the  four-point  derivative  approximation  formula.   Investigation 
showed  that  upon  using  the  coefficients  listed  by  Bickley,  the  resulting 
transfer  function  contained  an  additional  zero.   This  zero  was  of  the 
form  (St  +  l),  thus  it  would  cause  little  change  in  the  low-frequency 
response  of  the  circuit. 

While  Eq.  (^2)  gives  a  satisfactory  form  for  a  derivative  approx- 
imation, the  gain  requirement  dictated  by  the  6a, a  t   coefficient  re- 
quires further  consideration.   To  achieve  an  exact  derivative  in  the  low 
frequency  range  this  coefficient  must  equal  unity.   The   a  's   represent 
potentiometer  settings,  and  therefore  must  have  values  of  unity  or  less. 
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Since  the  a  a  's  represent  the  relative  gains  of  the  four  legs  of  the 
n  n      -^  to  to 

circuit,  the  gain  may  be  written  as 

(50) 


In  order  to  discriminate  against  60  cps,  as  with  the  other  circuits 
ssec 
Eq.  (50) 


discussed,  a  value  of   t  =  .005  sec.  was  found  convenient.   Hence  from 


-3  =  !^  (51) 

^   ^3 


which,  since 


^3^1 


gives 

OL^  >   900  (52) 

But  CC       was  defined  as 


"■s^^s 


and  therefore 

R^  >  900(r2  +  R^)  (53) 

Hence  it  follows  that  the  overall  system  gain  must  be  greater  than  900. 
A  gain  of  this  magnitude  makes  such  a  circuit  impractical  for  operations 
where  economy  dictates  simple,  unsophisticated  operational  amplifiers. 

The  theoretical  gain  curve  for  this  network  is  shown  in  Fig.  5-   As 
expected,  there  is  greater  accuracy  at  low  freq_uency  than  in  the  two- 
or  three-pole  cases  previously  discussed. 
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3.6.   The  Four-Pole,  Two-Leg  Derivative  Circuit. 

The  impracticability  of  the  previously  described  four-pole  system, 
arising  from  the  high  Inherent  amplification  requirement,  caused  con- 
sideration to  be  given  to  alternate  approaches.   The  first  of  these  is 
shown  in  Fig.  7-   If  for  this  circuit  we  set 

^1  -  ^I'-A 


'a   2r  +  B 


a,  = 


1   Sr^.R^ 


^3  =  "a^a^'a 


0L„    = 


2   2r^   +   R^ 


T,   = 


^2^2^2 


k        2v^   +  R^ 


Then  the  resulting  transfer  function  is  given  by 


e.  (s) 

m 


^A 


^2^2 


_(St^  +  1)(St2  +  1)  "  (St^  +  l)(STi^  +  !)__ 


(5^) 


This  tranfer  function  may  be  written  as 


eo(s) 
m 


.(St^  +1)(St2  +1)(St3  +  1)(St^  +  1] 


55) 


if  the  parameters  are  chosen  such  that  the  following  equations  are 
satisfied i 


^A  -  ^2^2  =  ° 


'^l'^2  -  "3^  =  ° 


(t^  +  T^)  -  (t^  +  T^)  ;^  0 


(56) 
(57) 
(58) 
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FIG.  7     FOUR  POLE,  TWO  LEG  DERIVATIVE  CIRCUIT 
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As  In  the  previous  derivative  approximations,  we  wish  to  choose 
the  time  constants  of  the  poles  such  that  60  cps  is  passed  without 
amplification.   This  requirement  provides  another  equation  relating  the 
T  's  ^  i.e., 

log  377  =  log  377t^  +  log  377T2  +  log  377T2  +  log  377^1^       (59) 

but  it  clearly  does  not  give  a  unique  solution.  The  optimum  approach 
appears  to  be  to  choose  two  of  the  x  values  and  determine  the  other 
two  from  these,  graphically.   Figure  8  shows  the  proceedure  adopted. 

The  graphics  of  Fig.  8  were  derived  from  the  high- and  low-frequency 
asympotes  of  the  transfer  function  given  in  Eq.  (55)-   The  procedure 
was  :   First  construct  the  four-equipole  envelope  with  o)  =  1  at  zero 
db  and  w  =  377  at  zero  db  as  the  predetennined  vertices  (the  lines 
OA  and  OB  have  slopes  of  +6db/octave  and   -l8db/octave  respectively) 
Second,  construct  the  double  pair  envelope.   This  was  done  by  choosing 
point  C  on  line   OA  at  a  position  corresponding  to  w  =  hO.      The  line 
CD  then  drawn  with  the  required  slope  of  -6db/octave,  hence  point  D 
was  determined.   Next,  through  point  C,  a  line  with  zero  slope  was 
drawn  and  through  point  D  a  line  with   -12db/octave   slope.   The  two 
undetermined  poles  of  the  four  unequal  pole  system  are  then  the  inter- 
cepts of  any  line  parallel  to  CD  on  the  lines  CE  and  ED. 

This  construction  method  is  completely  general  as  it  is  merely  a 
graphical  solution  to  Eqs.  (57)  and  (59)-   The  values  obtained  from  the 
example  given  in  Fig.  8  are: 

T^  =  .025        T^  =  .015 

T^  =  .005^1     T^  =  .009 

These  values  satisfy  all  the  requirements  of  Eqs.  (57)^  (58),  and 
(59)  and  yield  the  following  transfer  function  for  the  circuit  given 
in  Fig.  7- 
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Since   a   is  a  potentiometer  setting^   and  is  therefore  less  than 
unity,  for  the  correct  magnitude  of  the  derivative  at  low  frequencies, 
a   must  be  equal  to  or  larger  than  185.   This  gain  is  much  more  practical 
than  the  9OO  obtained  from  the  previous  four -pole  approximation  circuit. 

With  the  construction  method  outlined  in  Fig.  8,  however,  the  gain 
requirement  can  be  adjusted  since  it  is  in  effect  the  reciprocal  of 
[(t  +  T,  )  -  (t  +  T  )].   The  theoretical  frequency  response  of  this  net- 
work is  shown  in  Fig.  5  for  comparison  with  those  previously  discussed. 

All  circuits  have  been  examined  in  the  laboratory,  using  a  Donner 
Model  3500  Analog  Computer,  as  part  of  this  program.   Our  experiments 
show  that  the  output  of  each  circuit  for  various  control  type  inputs  and 
the  frequency  response  characteristics  are  in  agreement  with  our  theo- 
retical predictions.   The  results  of  the  frequency  response  tests  did 
not  vary  significantly  from  those  shown  in  Fig.  5"   The  response  of  the 
two-leg,  four-pole  derivative  device,  shown  in  Fig.  7^  to  a  range  of 
inputs  is  given  in  Figs.  9-13.   Results  of  the  other  circuits  showed  no 
marked  deviation  and  so  are  not  included.   Figures  9-13  show  the  input 
signals,  the  derivative  of  the  input  and  the  integrated  derivative  for 
comparison  with  the  input  signal. 

Following  the  method  of  synthesis  used  in  the  last  circuit,  several 
other  circuits  were  devised  using  various  combinations  of  R-C  time 
delays.   These  circuits,  along  with  their  respective  transfer  functions 
are  given  in  Appendix  A.   The  analysis  of  each  circuit  follows  the  for- 
mat previously  described. 

3.7.   An  Integral  Approximation  to  the  Derivative. 

In  Fig.  1^  we  show  a  comparison  between  a  Fourier  series,  a  Fejer 
sum  and  Lanczo's  CX     factor  curves  for  a  12  term  approximation  to  a 
rectangular  wave.   These  curves  indicate  that  the  Lanczo's  a  factor 
method  has  a  clear  advantage.   Although  we  shall  not  present  the  argument 
here,  we  point  out  that  the  Lanczo'  s  o;  method  is  essentially  the  same 
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approach  as  the  derivative  approach  to  which  we  have  previously  referred, 
It  seems  reasonable,  however,  to  consider  whether  or  not  there  may  be  a 
possible  analogy  with  the  Fejer  sum. 

If  we  refer  to  Fig.  15  then  we  see  that  the  derivative  is  approx- 
imately given  by 

(f^  -  f2)/2At  (61) 

Let  us  introduce  the  median  f  .   Then  the  area 

and  the  area 

Ag  ^  (fg  +  f^)  f^  (63) 


Thus 


or 


^  -  Ai  "  ('i  -  fs'  t"  (6^) 


.,2        2At       dt 
At 


(65) 


Now  the  area  A^   can  be  written  as  the  difference  of  two  integrals,  i.e., 
t-At  t-2At 


A     =    /          f(t)dt 

-J 

f(t)dt 

0 

0 

Likewise 

t 
p 

t-At 

A^  =    /  f(t)dt   - 

/ 

f(t)dt 

0 

0 

Therefore 

t 

t-At                       t-2At 

r                      r 

A^   -  A^=J  f(t)( 

it   -  2 

/       f(t)dt  +  /         f(t)dt 

0 

0                              0 

(67) 


(68) 
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Thus  it  follows  from  Eq.  (68)  that  the  derivative   df/dt  can  be 
expressed  in  terms  of  difference  in  integrals.   Thus 


i  -  (-)^ 


t-At  t-2At 


/  f(t)dt  -  2  /     f(t)dt  +  /      f(t)dt 


LO  0  0 


(69) 


We  examine  now  the  nature  of  this  approximation  to  a  derivative 
when  the  system  is  applied  to  a  signal  which  consists  of  two  components. 
A  smooth  signal  and  a  superposed  AC  signal.   We  see  immediately  that  if 
At  is  correctly  chosen  then  the  integrals  of  the  AC  components  are 
all  identical.   Thus,  since  the  algebraic  sum  of  the  gains  associated 
with  these  integrals  is  zero,  they  will  be  self -compensating  so  far  as 
AC   is  concerned.   We  observe,  too,  that  if  this  mathematical  expression 
is  to  be  transformed  into  an  electronic  circuit  then  the  input  signal 
must  be  trifurcated  and  the  three  arms  delayed  relative  to  each  other. 
At  takes  on  two  special  significances.   First,  it  is  equal  to  the  rel- 
ative time  delay  between  adjacent  arms  and  secondly  the  inverse  of  its 
square  is  the  overall  circuit  gain.   Thus,  as  before  the  system  does 
not  allow  complete  freedom  in  the  choice  of  At,  instead  At  must  be 
chosen  to  be  compatible  with  the  characteristics  of  the  amplifiers  which 
are  used.   The  integrating  amplifiers  act  as  very  definite  smoothing 
devices  for  cyclic  noise  and  the  delay  circuits  are  very  powerful  peak 
or  impulse  filters.   As  in  the  other  approximation!  circuits  the  smaller 
the  At  which  we  can  use  the  more  nearly  will  the  output  signal 
approach  the  true  derivative  level. 

3.8.   Integral  Difference  Derivative  Circuit. 

Figure  l6  shows  the  analog  network  used  to  implement  Eq.  (69).   The 
practicality  of  this  circuit  is  limited,  however,  because  of  the  tendency 
of  the  input  integrating  amplifier  to  saturate  for  any  input  of  constant 
polarity.   Analysis  of  the  transfer  function  for  the  network  given  in 
Fig.  16  shows  that  the  amplifiers  one  and  three  may  be  interchanged 
without  effecting  the  function.   Thus,  the  circuit  of  Fig.  17  gives  the 
same  result  as  that  of  Fig.  16,  but  prevents  saturation.   The  transfer 
function  of  this  circuit  is: 
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FIG.  16.    THREE  POLE  DERIVATIVE  CIRCUIT 
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e.  [s) 
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SRo^fo 


a  a 
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St^  +1   St^  +  1  ■  St^  +  1 


^2^2   ^  ^3^3 


(70) 


where 


T.  = 


R, 


a,  = 


r  +  R 
1    1 


T„  = 


r^  +  R^ 


R. 


o;^  = 


r  +  R 
2   2 


T„  = 


R. 


a,  = 


"•3*^3 


it  becomes 


e^s) 

e.  (sj 
m 


SRoCfo 


(St^  +  1)(St2  +  DCSt^  +  1} 


(71) 


if  the  parameters  are  chosen  such  that 


T  =  riT 

n 


^^1  -  ^2^2  ^  ^3^3  =  ° 


5^«1  -  ^^2^2  ^  ^^SS  =  ° 


6a^a^  -  3a2a2  +  Sa^a^  ;^  0 


(72) 
(73) 
(7M 
(75) 


The  solution  to  the  above  equations  for  A  a   yields  the  results 

n  n 

^2^2  =  2a^a^ 
^3«3  =  ^«1 


(76) 
(77) 


The  ratio  of  the  relative  gains  of  the  three  legs  is  l,-2,+l  and  these 
coefficients  are  those  given  by  Bickley  for  the  three  point  approximation 
to  a  second  derivative.   The  transfer  function,  Eq.  (7l)^  shows  the  network 

generates  a  second  derivative,  and  then  integrates  this  to  obtain  the  first 
derivative.   The  analysis  of  the  circuit  is  identical  to  that  discussed  pre- 
viously for  three  poles.   The  noteworth  point  of  this  approach  is  that  both 
the  first  and  second  derivative  are  generated  at  the  same  time. 
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IV.   CONCLUSIONS 


From  mathematical  approximations  of  the  type 


;,ni  (m) 
^r     1 


^Vo"  Vi"   ••••   VJ  '^  ^^2) 


mi    p:  '  o-'o   ri       p-'p 

A  systematic  family  of  derivative  circuits  can  be  derived.   These  circuits 
are  at  least  as  satisfactory  as  those  normally  adopted  and  in  most  cases 
superior  for  control  purposes.   Their  superiority  arises  from  the  fact 
that  they  incorporate  more  poles  than  is  customarily  achieved.   Thus  their 
"cut-off  frequency  line"  can  be  made  appreciably  steeper  than  the  con- 
ventional -6db/octave.   The  slope  of  this  line  is   -6(n  -  l)db/octave 
where  n  =  pole  number.   Associated  with  this  there  is  an  increase  of 
accuracy  in  derivative  value  for  low-frequency  signals. 
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